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ABSTRACT 

We detect emission from multiple low-excitation ro- vibrational transitions of 
OH from the two Herbig Ae stars AB Aurigae and MWC 758 in the 3.0 — 3.7/xm 
wavelength range (L-band), using the NIRSPEC instrument on Keck II. The 
inner radius for the emitting region in both stars is close to 1 AU. We compare an 
optically thin LTE model and a thin-wedge fluorescence model, finding rotational 
temperatures of 650 — 800 K and OH abundances of 10*^^ — lO''^ molecules for the 
two stars. Comparisons with current chemical models support the fluorescence 
excitation model for AB Aurigae and possibly MWC 758, but further observations 
and detailed modeling are necessary to improve constraints on OH emission in 
different disk environments. 

Subject headings: stars: pre-main sequence — planetary systems:formation and 
protoplanetary disks — infrared: stars — molecular processes — circumstellar 
matter 



1. Introduction 

We report the first detection of warm hydroxyl (OH) emission from the circumstellar 
disks of two young intermediate-mass stars. AB Aurigae (HD 31293) has a mass of 2.4 M© 
and an age of 1-3 Myr, while MWC 758 (HD 36112) has a mass of 2.0 Mq and an age of 2-5 



iNASA Goddard Space Flight Center, Code 693, Greenbelt, MD 20771, USA 

•^California Institute of Technology, Division of Geological and Planetary Sciences, MS 150-21, Pasadena, 
CA 91125, USA 

■^Catholic University of America, Department of Physics, Washington, DC 20064, USA 
^Corresponding Email: Avi.Mandell@nasa.gov 
^NPP Postdoctoral Fellow 



- 2 - 



Myr (Ivan den Ancker et al.lll998l ). Each disk is seen at low inclination, implying that our 
observations probe the surface layer and warm inner region {R < 20 AU) of the circumstellar 
material. 

The detection of OH provides a new window into the gas chemistry of these disks. OH 
has a central role in many chemical reactions involving the formation and destruction of H2, 
CO, and H2O, and the detection of warm hydroxyl opens up an excellent new diagnostic of 
the thermal and radiative environment in the planet formation zone of circumstellar disks. 
We achieve very high S/N results using improved data reduction algorithms, allowing us to 
sensitively test different excitation models. These results, combined with similar data for 
an expanding suite of molecular tracers such as CO and H2O, will significantly improve our 
ability to constrain the disk chemistry driven by intense radiation from the central star. 



Observations and Data Reduction 



Spectra were acquired at high resolving power (X/A X ~ 27,000) using t he Near-IR 
Spectrograph (NIRSPEC) on the 10-m Keck H telescope (IMcLean et al.lll998l ). P-branch 
transitions {J" = 2.5 to J" = 17.5) in the OH (1-0) band were sampled on observing runs 
in April, October, and December 2006. Each science target was observed on all three runs; 
the comparison star HR1620 was also observed as close in time to the two science targets 
as possible (on the December run a second comparison star (HR2714) was observed as an 
additional check on the reduction methods). The telescope was nodded 12 arcsec in an 
ABBA sequence, with 60 seconds integration per beam. 

Spectra were extracted and p rocessed for each echelle order in eac h ABBA set using 
custom data reduction algorithms (IMandelll 120071 : IVillanueva et al.ll2008l ): in this Letter we 
provide a summary of these procedures, with further details to be presented in a follow-up 
paper. The terrestrial spect ral transraittarice was synthesized for each ABBA set using the 
GENLN2 atmospheric code (lEdwardslll992l : iHewagama et al.ll2003l ) and subtracted from the 
stellar spectrum; the difference removed the stellar continuum and compensated for changing 
airmass and telluric water vapor. We co-added the residuals for all sets of each star, and 
the mean residuals of the science and comparison stars were then differenced, removing 
remnant fringes and other instrumental artifacts. This process achieved photon-noise-limited 
results corresponding to S/N ~ 2000 on the original stellar continuum (see Figured]). Flux 
calibration was performed by comparing the observed continuum flux near each transition to 
the predict ed flux based on th e L-ba nd stellar magnitu de and temperature. Magnitudes were 
taken from lHillenbrand et all Jl992h (AB Aurigae) and lMalfait et all Jl998l ) (MWC 758); for 
both stars, we used temperatures derived from fits to the NIR spectral energy distribution 
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bv lNatta et all fl200lh . 

The measured shapes of detected hnes are asymmetric to varying degrees because they 
remain convolved with the atmospheric transmittance function; other unknown factors may 
also contribute. The strongest OH doublet detected (P4.5 1+,1-) is not severely affected by 
local terrestrial lines and the line shapes are adequately reproduced by a Gaussian emission 
kernel convolved with an instrumental line shape function (Gaussian, with parameters deter- 
mined from measured atmospheric features). We solved for the "best-fit" velocity width and 
radial velocity offset for the emission kernel, and assumed these values for other OH lines. 
We then corrected the line areas measured for other OH doublets by fitting the data and 
measured errors to the scaled Gaussian line shape model convolved with the atmospheric 
transmittance in the local spectral region and the instrumental line shape; the uncertainty 
for each line was calculated using the standard errors for the derived model parameters of 
the fit. 



3. Results 

For both stars, either one or both A-doublet components were detected for all available 
OH V = 1 — ro-vibrational transitions ranging from J" = 2.5 up to J" = 9.5 and excitation 
energies up to 1600 cm~^. Scaled line fluxes for both stars are shown in Figure [21 individual 
fluxes and uncertainties for each line will be tabulated in a follow-up paper. Lines were 
detected in all three sets of observations, and for each star the detected lines were found to 
have a common radial velocity shift (within 3 km s~^) relative to the barycenter. Radial 
velocities derived from fits to the strongest and clearest doublet (P4.5 !+,!-) are 16.9 ± 1.1 
and 15.7 ± 0.7 for AB Aurigae and MWC 758 respectivel y; these values co mpare well with 



radial velocities derived from sub-mm CO observations (iDent et al.ll2005l ). At NIRSPEC 
resolution (13 km s~^) the line shapes of the strongest doublet (P4.5) are only marginally 
resolved; the FWHM for the P4.5 components (iteratively deconvolved from the instrumental 
profile as described in |J2]) is 16.2 ± 2.8 km s~^ and 22.1 ± 1.5 km s"^ for AB Aurigae and 
MWC 758, respectively. 

The measured line fluxes can be related to an excitation temperature Tex and total OH 
abundance Nt^t using models for the excitation physics and radiative transfer. The simplest 
model assumes a coUisionally-dominated optically thin medium in which the vibrational and 
rotational distributions for OH are described by a single temperature (LTE). However, if the 
emitting OH is located in a region of the disk atmosphere where the timescale for collisional 
excitation is longer than the timescale for radiative excitation of ground-state OH, rotational 
state distributions in the v = 1 vibrational level is instead governed by both the ground state 
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populations and the incident radiation field. We consider a simple collisional excitation model 
in §3.11 and a radiative excitation model in §3.2( non-LTE collisional effects are deferred for 
future consideration. 



3.1. Collisional Excitation Model 

When collisions dominate the excitation process, the emission intensity of ro- vibrational 
transitions varies lin early with the number of molecules and logarithmically with the exci- 



tation temperature (jHerzbergJ Il950l ) . A rotation diagram, graphing the mean flux for each 



doublet pair (normalized by uAw') against the energy of the excited level (J', v = 1), is shown 



in Fi gure [21 transition parameters are taken from the HITRAN database (IRothman et al. 



20051 ). A linear fit to the rotation diagram for each star yields the excitation temperature 
(Tea;), and we use the E' = intercept from each fit multiplied by the total partition func- 
tion for the derived tempera tures (from HITRAN) and scaling by the distance to each star 



(Ivan den Ancker et al.l 119981 ) to find the total number of OH molecules. Derived values are 



listed in Table [H 

An inner radius (/?«„) for the OH emitting region can be found if the emission lines are 
broadened primarily by the projected radial velocities of material in Keplerian orbits around 
the central star. We assume a maximum orbital velocity of 2xHWHM of the P4.5 lines, after 
removing the contribution from thermal broadening (~0.8 km s^^ for OH at 700K). However, 
the derived inner radius varies with the local disk inclination (i) as sin^(i), and published 
inclinations for circumstellar disks from NIR interferometry can vary substantially. For 
MWC 758, values from Eisner et al. (2004, hereafter E04) and Isella et al. (2006, hereafter 
106) are similar, but for AB Aurigae values of sin(i) derived by E04 and 106 differ by more 
than a factor of 2. Values for Rin calculated using both sets of inclinations are given in Table 
[U an outer radius for the emitting region is not constrained by the line shapes due to the 
low disk inclinations and limited spectral resolution. 



3.2. Radiative Excitation Model 

Fluorescent OH emission can occur in diffuse astrophysical environments that are strongly 
irradiated, but the excitation details depend critically on the local UV and IR radiation field. 
OH can be excited to the v = 1 vibrational level (in the ground electronic state) both through 
ro-vibrational transitions in the near-IR as well as electronic transitions in the UV followed 
by cascades to the ground electronic state. OH fluorescence has been studied extensively in 
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comets (ISchleicher fc A'Hearnlll982l : IWeaver fc Mummalll984l : ISchleicher fc A'Hearru 1988 ) 



and UV fluorescence of CO lias recently been explored in disks (IBrittain et al 



20031 . l2007f ). 



We assess the effectiveness of OH fluorescence through a fit to a simple thin wedge 
model, with free parameters being the scaled thickness of the wedge (H/R), the column 
density in the absorbing (radial) direction at the inner boundary {Ncoi,r), and the lower-state 
rotational temperature (Tiow) (Figure [3]). Integrating over the wedge geometry conserves the 
total number of available photons with distance from the star, and the use of an absorbing 
column density removes the need for a radial density profile. For each Ncoi,r, an optical depth 
and equivalent width is calculated for each transition and multiplied by the radiative flux at 
the minimum orbital distance {Rin, as defined in §3.11) to find the total absorbed flux. Model 
parameters {Tiow and H/ R) are constrained by comparing predicted and measured fluxes for 
each transition; the value for the fit is used to evaluate the effects of variations in optical 
depth for different transitions. Multiplying Ncoi,r by H at the inner boundary gives the total 
number of molecules Ntot- 

The fluorescent emission per molecule is calculated by summing the pump rates for 
transitions from all levels that contribute to state J' using the Einstein B coefficients and 
then calculating the fraction of elec trons that decay in to state J" using branching ratios 
based on the Einstein A coefficients (jPiSanti et al. 2006 ). For UV fluorescence, we consider 
pumping into the three lowest vibrational levels in the first electronic state (A^S) and decay 
into v" = 1 in the ground state (X^II) using rotational transition pa rameters derived from 
publi shed vibrational band frequencies and transition probabilities (ISchleicher fc A'Hearn 
19881 ). UV fluxes in the appropriate spectral region are taken from lUE spectra for both AB 
Aurigae and MWC 758 (jRodrfguez-Pascual et al.lll999l ): our approach does not compensate 
for flux variat ions on small spectral scales that cause UV fluorescence variations in comets 
( 1Swingslll94ll ) but is adequate for this study. For IR fluorescence, we calculate pumping and 
decay rates using parameters from HITRAN. We adopt a basic 2-component model for the 
IR radiation field, comprise d of thermal emiss ion from both the star and an annular wall of 
"hot dust" surrounding it (iNatta et al.ll200ll ): comparisons of our measured L-band fluxes 
for AB Aurigae and MWC 758 with those of standard stars of the same spectral type suggest 
the "wall" emits approximately 10 times more L-band flux than the star alone. 

The location of the NIR-emitting dust region for both stars is uncertain. SEDs of 
HAe/Be stars are fit best by dust with temperatures of ~1500K (t he approximate su blima- 
tion temperature for "astronomical silicate" dust) for the inner rim (INatta et al.ll200ll ). If we 
assume the dust has an equihbrium temperature determined solely by the stellar radiation 
field, the dust rim would be located at a radius of 0.48 AU for AB Aurigae and 0.33 AU 



for MWC 758 (using stellar luminosities and stellar radii tabulated by Ivan den Ancker et al. 
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( 119981 ) )■ These values are consistent with the analysis of NIR interferometric data by 106; 
however, earlier work by E04 found significantly smaller inner dust radii for the two stars. 
For completeness, we ran models using inner dust radii and disk inclinations (important for 
defining Rin) from both studies; model predictions are tabulated as a range of values for 
each star in Table [H 

Predictions of our thin wedge model for various parameters (slab thickness (H/R), 
lower-state temperature (Tiow), and x^); iterated over different values of the absorbing OH 
column density, are shown for AB Aurigae in Figure [31 A Ax^ of less than 9 is predicted 
by either optically thin or optically thick absorbing column; however, optically thin models 
have H/R > 10~^ while optically thick models have H/R < 10~^. Predictions for Tiow and 
from the optically thin models (assuming a slab thickness limit of H/R > 0.01) are listed 
in Table [H We do not consider the optically thick models realistic as they would require 
\og{Ntot) > 51 and imply an enormous gas:dust ratio and hydrogen abundance, but this 
scenario cannot be ruled out with the current data set. 



Discussion 



The primary pathway for OH production considered in models of irradiated environ- 
ments inv olves the recombination of O and H atoms after destr uction of H2 and CO by UV 
radiation (IStoerzer fc Hollenbachlll998l ) or collisional processes ( iKamp fc Dullemondll2004l ) . 
Excited OH can also be forme d directly from t h e photolysis of H2O , producing emissio n 
known as "prompt emission" (ICarringtonl ll964l : iMumma et al.l ll97ll : iBonev et al.l 120061 ) . 
However, in comets the OH responsible for prompt emission is not characterized by a single 
rotational temperature, so it is most likely not responsible for the emission reported in this 
paper. Understanding the role of these OH formation and excitation processes and where 
they occur in circumstellar disks has the potential to provide a number of constraints on 
disk structure and chemistry. 

Th estimates for our derived inner radii using the most recent interferometry results 
{Rin ~ 1 AU) and our derived temperatures (T ~ 700K) suggest the detected OH emission 
originates in the surface layers of the disk. The line shapes for OH are equivalent within 
error to those found for low-J CO v = 1 — ro-vibrational transitions (IBlake fc Boogert 



20041 ). suggesting a common inner radius. 



Our results for rotational temperature and OH abundance can also be co mpared with 



recent chemical models of circumstellar disks. De-coupled gas-dust models by iKamp et al. 



( I2OO5I ) predict high disk temperatures (T > 500K) just above the optically thick region of 
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the disk midplane owing to the confluence of radiative and viscous lieating, resulting in tlie 
production of a population of warm OH in the disk atmosphere; below the optically thick 
dust boundary OH production ceases. For AB Aurigae, this model predicts vertical column 
densities of N(.oi,v ~ 10^^ mol cm~^ for the warm OH, extend ing from the inner disk out 



to 15 AU (I. Kamp, private communication). For a 1 star, iMarkwick et al.l (120021 ) also 
predict OH column densities of 10^^ — 10^^ mol cm~^ at 10 AU; however, the location and 
temperature of the OH produced in these models is uncertain. If we assume an outer radius 
of 15 AU and a constant vertical column density, results from our fluorescence model imply 
NcoLv ~ 10^^ mol cm"^ for AB Aurigae and Ncoi,v ~ 10^^ mol cm"^ for MWC 758. Results 
from our LTE model would be two orders of magnitude larger than the model predictions 
for AB Aurigae, but are less divergent for MWC 758. 



In concert with this study, ISalyk et al.l (120081 ) reported a detection of OH and H2O in 



two T Tauri stars, with inner radii and temperatures similar to those found here. They use 
optically thick CO lines to find an emitting area (0.2 — 0.4 AU^); converting their column 
densities to total molecules gives \og{Ntot) ~ 43 based on an LTE excitation model. These 
abundances are consistent with ours; additional analysis of the two samples will be presented 
in future work. Sampling the same spectral region, we do not detect H2O emission lines 
in either AB Aurigae or MWC 758. Using the S-cx uncertainty for the spectral region as 
an upper limit for a line detection and applying the LTE analysis from above, we derive 
preliminary upper limits on the abundance of emitting H2O that are lower than our derived 
OH abundance (see Tabled]). HAe stars have high UV fluxes that may dissociate H2O in the 
disk surface, but a detailed model of the excitation and dissociation physics of water under 
these conditions has not yet been constructed. 
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Table 1. Quantities Derived in This Work 



Model AB Aurigae MWC 758 

Inner OH Radius'' E04 O.U+q I^ AU 1.29l°;44 AU 

106 l-05t[;-jg AU 1.62t°-Jf AU 

LTE Tex 651K±41K 746K±39K 

OH: log{Ntot) 44.20±0.28 43.73±0.28 

H2O: log{Ntot) < 44.4 < 43.7 

Fluorescence T„^o 690K-725K 751K-808K 

OH: \og{Ntot) 42.0-42.9 43.5-44.5 



Ca lculated using disk inc linations from lEisner et al.l 
l|2004l) (E04') and llsella et al. 
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Fig. 1. — Examples of line detections for both stars. The upper section of each frame shows 
the stellar spectrum with terrestrial lines, as measured; a synthetic atmospheric spectrum is 
overlaid for comparison. At the bottom of each frame are individual residuals for the science 
target and the comparison star after removal of a model atmosphere. The middle residuals 
are the difference of the science and comparison star residuals bounded by 3-cr photon-noise 
uncertainty limits. A S/N ~ 2000 was achieved in all spectral regions without severe telluric 
absorption. 
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Fig. 2. — Rotation diagrams of detected transitions for AB Aurigae and MWC 758, with 
both doublet components averaged to minimize uncertainties. We detected all transitions 
with J" < 9.5 having transmittance above the noise limit; for two transitions (P3.5 (2), 
P6.5 (1)), only one component was measurable due to the atmospheric transmittance and 
therefore only the lower limit is secure. The P9.5 (2+, 2-) transition was not included in the 
fit for AB Aurigae; the inconsistently high flux may be due to a blend or an instrumental 
effect. 
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Fig. 3. — Predictions of a thin wedge fluorescence model for AB Aurigae. For a specific 
absorbing column density, a lower-state rotational temperature (Tiow) and wedge thickness 
{H/ R) are derived to match the observed component-averaged fluxes for each transition; the 
Ax^ value (compared to the minimum value of 1.52) and standard confidence intervals for 
the fit are also shown. 



